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ABSTRACT 

Context. NGC 6791 is a unique stellar system among Galactic open clusters, which is at the same time one of the oldest open clusters 
and the most metal rich. It is located inside the solar circle, harbors a large population of binary stars, and possibly experienced 
prolonged star formation. The combination of all these properties is puzzling and poses the intriguing question of its origin. 
Aims. One possible scenario is that the cluster formed close to the Galactic Center and later migrated outward to its current location. 
In this work we study the cluster's orbit and investigate the possible migration processes that may have displaced NGC 6791 to its 
present-day position, under the assumption that it actually formed in the inner disk. 

Methods. To this aim we performed integrations of NGC 6791 's orbit in a potential consistent with the main Milky Way parameters. 
In addition to analytical expressions for halo, bulge and disk, we also consider the effect of bar and spiral arm perturbations, which are 
expected to be very important for the disk dynamical evolution, especially inside the solar circle. Starting from state-of-the art initial 
conditions for NGC 6791, we calculated 1000 orbits back in time for about 1 Gyr turning different non-axisymmetric components 
of the global potential on and off. We then compared statistical estimates of the cluster's recent orbital parameters with the orbital 
parameters of 10^ test-particles originating close to the Galactic Center (with initial galocentric radii in the range of 3-5 kpc) and 
undergoing radial migration during 8 Gyr of forward integration. 

Results. We find that a model that incorporates a strong bar and spiral arm perturbations can indeed be responsible for the migration of 
NGC 6791 from the inner disk (galocentric radii of 3-5 kpc) to its present-day location. Such a model can provide orbital parameters 
that are close enough to the observed ones. However, the probability of this scenario as it results from our investigations is very low. 

Key words. Galaxy: kinematics and dynamics - open clusters and associations: general - open clusters and associations: individual: 
NGC 6791 - Galaxy: structure 

1 ■ Introduction abundance of iron and other elements much larger than observa- 
tional uncertainties (Milone et al. 2009) . 

Star clusters -globular and open -have long been recognized as The rich, intermediate-age star clusters of the Magellanic 

ideal benchmarks for stellar evolution theories primarily because Clouds (MCs) do not seem to be an exception ( |Milone et al. | 

they are made of a large number of stars that share the same age, ^009), although m this case signatures of prolonged star forma- 

^ ■ distance, and metallicity. This allows comparisons with stellar tionaremorecommonthanevidencesof discrete episodes of star 

^- models in a statistical sense. The only exception has routinely formation. Indeed, m the CMD of MCs clusters the MS is much 

been oj Gen, which has been known for quite some time to dis- broader than expected by photometric errors and bmary stars, 

play a significant spread in color in its red g iant branch (RGB) but does not separate into discrete sequences unlike m Galactic 
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and inhomogeneities in chemical elements (JPiotto et al.l l2QQ5l giobulars. 

IVillanovaetaDl2QlQl and references therein). As of now, the only MW open cluster where evidence of 

prolonged star formation has been found is NGC 6791 . The up- 

However, in the past decades, signatures of multiple stellar per MS and turn-off region of this cluster are much wider than 

populations - with different age and metallicity - have been de- expected from photometric errors and binary stars alone. By 

tected in virtually all Milky Way (MW) ^lobul ars by means of separating the stars lying inside the cluster core and outside it, 

very accurate photometry and spectroscopy (see |Piottoll2QQ9l for iTwarog et al.l (l201ll) found that inner stars occupy the red edge 

a review on the subject). Typical photometric signatures are the of the upper MS, while outer stars occupy the blue edge. The 

presence of multiple distinct main sequences (MS) or RGBs in most viable explanation for this occurrence is an age spread of 

the same cluster color magnitude diagram (CMD), while spec- ~ 1 Gyr between the inner and the outer regions of the cluster, 

troscopic signatures are normally in the form of spreads in the in the sense that inner stars formed before outer stars. Star for- 
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mation therefore started probably in the center of the cluster and 
continued outward during quite a long time. 

NGC 6791 is unique in the MW because it is one of the 
most massive old open clusters with an age of ~8Gyr, and 
is extremely metal -rich ([Fe/H] ~ +0.40. ICarraro et al.l |2006[ 
lQrigliaetal]l2006h . No other open cluster in the MW is that 
rich in metals. This combination of age and metallicity is unique 
in the whole population of Galactic open cluster s. While the 
RGB of NGC 6791 is broad in color (|Janeslll984 . no statisti- 
cally significant evidence for a metallicity spread within the clus- 
ter -which would support a self-enrichment scenario -has been 
reported so far. Within observational errors, stars in NGC 6791 
show t he same metal abundance. As mentioned by lTwarog et al.l 
(1201 ll) . since the CMD-hased age spread estimated for popu- 
lous clusters is usually less than the 0.7 Gyr found for NGC 419 
by\kubele et al {2010), the CMD of NGC 6791 may have resem- 
bled that of the Magellanic Cloud clusters 7 Gyr ago. 

The cluster's current location is at 4kpc from the Sun, close 
to the solar circle, at 8 kpc from the Galactic C enter (GC) and 
about 0.8 kpc above the plane (derived from Brogaar d et al.l 
120111 The metallicity of this cluster is at odds with the ra- 
dial abundance grad ient defined by the other MW open clusters 
(ICarraro et al.ll2006h . 

This wealth of unique properties clearly poses the basic 
question of the origin of this system. So far, to answer to this 
question only arguments based on the observed chemical charac- 
teristics of the cluster have been used. If this system is so metal 
rich, it must have formed very close to the Galactic bulge and 
then been deposited where it is now by some migration process. 
Indeed, only in the inner regions of the MW can star forma- 
tion and metal enrichment be fast and strong enough to generate 
a stellar system like NGC 6791. 

Alternatively, one can argue that NGC 6791 is the left- 
over of a dwarf elliptical that was disrupted by the MW some 
time ago and was then engulfed by the Galaxy, losing most 
of its mass by tidal interaction. This seems to be supported 
by the many e xtreme blue horizontal branch stars discovered 
in the cluster (iKaluznv & Rucinskil 1 19951 : iLiebert et all 1 19941: 
iBuson et aP l2006l) and by the recent measure of a significant 
UV-upturn in the sp ectral energy distribution of the cluster 
(iBuzzoni et al.ll2012b . Such a UV-upturn is only found in mod- 
erately UV-enhanced ellipticals and has the same strength as in, 
e.g., NGC 4552 and NGC 4486, which emit ~ 2% of their light 
at a shorter wavelength than 2500 A. This feature does not have 
any counterpart in the population of Galactic open clusters. 

The study of the NGC 6791 Galactic orbit did not reveal 
much to clarify the scenario so far, since depending on the qual- 
ity of present-day kinematic data, simple backward calculation 
could produce diff'erent orbits - a regular or bit wi th di fferent ec- 
centric ities e ~ 0.5 or ~ 0.3 bv .Bedin et aP (l2006l) . and lWu et aD 
(l2009l) . respectiv ely, both reach i ng gal ocentric radii of ~ 10 kpc, 
and an orbit by ICarraro et al.l (l2006l) extending to more than 
20 kpc in the Galactic disk. 

In this paper we present a numerical study of the NGC 6791 
Galactic orbit that includes the combined effect of the Galactic 
bar and spiral arms into an analytical MW model. At the same 
time, we aim at investigating in detail the scenario in which 
NGC 6791 formed close to the GC, where star formation might 
have been rapid and efficient enough to generate it. This scenario 
requires some radial migration process to move the cluster' s or- 
bit outward in the Galactic disk. 

Several different dynamical mechanisms triggering radial 
mixing are discussed in the recent literature: 



1 . the resonant sc attering by transient spiral st ructure was firstly 
introduced by ISellwood & B inney (zOOZ). Stars can move 
inward and outward at the corotation radius while preserv- 
ing their circular orbit. During the evolution of the disk the 
angular velocity (the corotation radius) changes and the spi- 
rals can induce radial migration within the whole disk. The 
importance of the transient spiral structure radial migration 
phenomenon was confirmed by Roskar et al. (2008 ) (see also 
Roskaretal. 2011). 

2. iMinchev & Famaeyl (1201 Ol) described another diffusion 
mechanism resulting from the resonance overlap of multi- 
ple rotating patterns, the Galactic bar, and spiral arms. Using 
test-particle simulations, these authors found that the over- 
lap mechanism can be up to an order of magnitude more ef- 
fective than the one due to transient spiral structure (for the 
model consistent with the MW it takes ~ 3 Gyr to produce 
the same radia l mixing as in ~ 9 Gyr by the transients, see 
also iMinchev & Famaeyl 120111) . This migration mechanism 
is a strong function of the strength of MW bar and spiral 
arms. However, this dependency is non-linear. The charac- 
teristics o f the resonance overlap mechanism were later con- 
firmed by IMinchev et al.l (1201 ll) via a fully self-consistent, 
tree-SPH N-body simulations. Among other results, these 
authors also found that the presence of a gaseous component 
makes the migration even more efficient. 

3. Small satellites ha ve been found to effectively mix the outer 
disks of galaxies. lOuillen et al.l (12009 ) used disk-like dis- 
tributed test-particles perturbed by a low-mass satellite (of 
mass a few times 10^ Mq). These authors find that the disk 
stars from the outer disk can m ove to low e ccentricity or- 
bits inside their birth radii. Bir d et al.l (1201 2l) compared ra- 
dial migration in several fully self-consistent N-body sim- 
ulations both with and without interaction with satellites. 
Among other conclusions, these authors showed that the ra- 
dial migration induced by satellites is a distinct mechanism 
that does not happen in isolated systems. 

During the past decades it has become widely accepted 
that the MW has a central bar. Evidence for this comes 
from observatio ns of t he inner Galaxy (e.g ., NIR photometry - 
Blitz & Spergel I1991L Dw eketal. 1995 1 iBinney et al.1 119971: 
gas kinematics -Englmaie r & GerhardI U999; and microlens- 
ing - Rattenburv et al. 2007 ) and dynamical effec ts near the solar 
circle (Dehnen 2000; Minch ev et al]l2007ll2010l) . It is, therefore, 
very important to consider the MW's bar in orbital integration 
studies, especially i n the inner disk, but a lso nea r the Sun. 

While for the ISellwood & Binneyl (l2002l) mechanism to 
work short-lived transient spirals are required, in barred galaxies, 
such as the MW, spiral s are more likely coupled with the bar, and 
are th us longer-lived (iBinney & Tremainel 12008; Ouil len et al.l 
I2OIOI) . Therefore, in this paper we will concentrate on the mi- 
gration induced by spiral-bar interaction, rather than transient 
spirals. The migration induced by orbiting satellites will also be 
unimportant for the inner disk. 

In Sect. [2] we describe the adopted analytical Galactic model 
and introduce spiral arm and bar perturbation. The initial condi- 
tions and the results of the backward orbit integration are pre- 
sented in Sect. O We concentrate on radial migration in Sect. |4l 
and, finally, discuss the results of this work in Sect. |5] 

2. Galactic models and orbit integrator 

To compute any Galactic orbit, it is necessary to adopt a model 
for the MW gravitational potential. The potential should repro- 
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duce the observed mass density of the MW. While the general 
picture of the today's MW structure is more or less well es- 
tablished and supported by observations, a m ore detailed de- 
scripti on is still under lively discus sion, see e.g. JSchonrich et afl 
(l2Qld) or lCoskunoglu et aP (l2Qllh for the rec ent analysis of the 
velocity of the local standard of rest (LSR), or lReid et all (l2009l) 
for the up-to-date description of MW rotation curve. These ques- 
tions naturally bring difficultly quantifiable uncertainties for the 
construction of any MW model. 

We chose to use a relatively simple model composed of sev- 
eral components. The axisymmetric time-independent part was 
modeled by a bulge, disk and halo. Additionally, we included 
a non-axisymmetric component -bar and spiral arms -that can 
act separately or together. For each component, an analytic ex- 
pression for the density is provided, from which the correspond- 
ing potential and force can be derived. 

The model is described in the following subsection. The 
main idea of our study is to qualitatively estimate the viability 
of the migration induced by the bar and spiral arms resonance 
overlap mechanism relative to the migration in simpler poten- 
tials, either purely axisymmetric or containing just one type of 
non-axisymmetric perturbation. 



2.1. Axisymmetric model 

For the time-independent axisymmetric compo nent of the 



Galact ic potential we used the model introduced bylFlynn et al.l 
(Il996h that was also used by [Gardner & Flvnnl (l2QlQh . V\ This 
model has three components - bulge, disk, and dark halo. The 
bulge is modele d as a superposition of two Plummer spheres 
(|Plummejl91lh . The dis k is modeled as a superpos ition of three 
Miyamoto — Nagai disks (iMiyamoto & Nagail 19751) . Finally, the 
dark halo is modeled as a spherical logarithmic potential. For 
clarity, we give the equation for each component of the axisym- 
metric potential: 

(1) 

(2) 



GMc, GMc, 



V^^^ 4 



2 I r2 



yZ _|_ y. 



C2 



0D„ 



-GMd 



^{7?2 + [a,+ Vfe' + ^')f} 



n= 1,2,3, 



0H = \vl\n{? + ?^. 



(3) 



(4) 



where r is the galocentric radius (r^ = 7?^ -h z^), G is the gravi- 
tational constant, M^ and Mq^, ^Ci and rc2, are the masses and 
scale radii of the bulge components, respectively. Md„ and a„ are 
masses and scale-lengths of the three disk components {n goes 
from 1 to 3). The parameter b is related to the scale-height of 
the disk and is the same for all three components of the disk. 
The halo has core radius ro and Vh is the asymptotic circular 
velocity at large radii (relative to ro). 



^ We also tested the model of lAllen & SantillanI (11991 '). slightly mod- 
ified to follow the more recent observations rotation curve as described 
below. The nature of the general results presented in Sects. [3] and |4] 
stays similar as in the case when the axisymmetric component is mod- 
eled as d escribed in Sect. 12.11 In its simplicity the model introduced by 
I Allen & S antillan (1991) does not follow the most updated knowledge 
of the MW mass distribution. However, it has been recently and exten- 
sively us ed by several teams dealing with orbits of MW clusters (e.g., 
lAllen et al. 2006, 2008; Bellini et al. 2010; Wu et al. 2009) or with the 
solar neighborhood kinematics (Antoia et al.ii2009i.i20 ID . 



Table 1. Parameters of the axisymmetrical component of today's 
MW model. 



Component Parameter 



Bulge 



Disk 



Value 



Halo 



Mci 


U.J-iU^"MO 


^Ci 


2.7 kpc 


Mc, 


\.6-WMq 


rc2 


0.42 kpc 


Md, 


8.9-1010Mo 


ai 


5.0 kpc 


Md2 


-6.9-1010Mo 


^2 


15.8 kpc 


Md3 


2.8-1010Mo 


(23 


33.0 kpc 


b 


0.3 kpc 


Vh 


225kms-i 


^0 


8.4 kpc 



iFlynn et aP (Il996l) chose the parameters of their model to 
reproduce the current MW observations at that time. We mod- 
ified the values of some parameters to achieve a model consis- 
tent with more recent MW observations. As for the recent ob- 
serv ations of the rotat ion curve at the solar galocentric distance 
Rq, iReid et al.1 (l2Q09h estimated Rq = (8.4 + 0.6) kpc and the 
solar circular rotation speed 0© = (254 ± 16)kms"^ To de- 
rive these values, they used the solar motion with respect to the 
LSR, v© LSR' ^s determined by iDehnen & Binneyl (1 1998b from 
the Hipparco s stars. The sanie Hipp arcos data set was later re- 
examined by lSchonrichet al.l (l2010h . who obtained slightly dif- 
ferent values for the vqlsr components. Namely, they found 
the velocity component in the direc tion of the Galact i c rota - 
tion to be about 7 km s"^ higher than JDehnen & Binneyl (1 19981) . 
Cons idering the s e mor e recent results, we changed the orig- 
inal Flynn et all (Il996l) values for the halo circular velocity 
Vh and the disk parameters (masses and scale radii) to obtain 
00 = 243 km s"^ at Rq = 8.4 kpc. These values are consis- 
tent with the &(7) derived from the apparent proper motion of 
SgrA* (.Reid & Brunthaie3l2QQ4 ) using Rq = 8.4 kpc and the 
ISchonrich et al.l (l2Q10h solar motion. 

The superposition of three Miyamoto — Nagai disks builds 
up a disk with exponential falloff' of the radial surface densit y 
on a wide range of radii. Similar as in [Gardner & Flynnl(l2QlQl) . 
we chose the disk parameters to obtain the scale-length of the 
exponential falloff' of the disk consistent with recent observa- 
tions. Our disk gives a radial surface density profile approxi- 
mately corresponding to the exponential falloff' with the scale- 
length of 2.6 kpc fo r radii of 5-18 kp c, which is consistent with 
recent observations ("Ju ric et aDl2Q08h . 

Values for all parameters of the axisymmetrical component 
are given in Table [T] 

The above described axisymmetric potential was also used 
in the models including the non-axisymmetric components. 
However, when we included the bar perturbation, the mass Mc2 
of one of the bulge components was reduced and the bar re- 
placed most of the bulge mass. We kept the axisymmetric com- 
ponents unaltered when we turn the spiral arms perturbation on 
(see Subsection 12. 31) . 



2.2. Galactic bar 

For the bar model we adopted the Ferrers potential (see para- 
graph 2.5.3 (a) of Binney & Tremaine 2008, and we chos e 
n = 2) of an inhomogeneous triaxial ellipsoid (|Pfennigej|l984b . 
The bar model is defined by the bar length, axes ratios, mass. 
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Table 2. Parameters of the non-axisymmetrical components in 
today's MW model (indicated as MWl in Sect. ^ with their 
references. 



[FR,tot(x.y)-<FR,tot(R)>]/<FR,tot(R)> 



F<^(x,y)/<FR,,ot(R)> 
•0.3 -0.2 -0.1 



Parameter 



Value 



Reference 



Bar mass 

Bulge mass Mc2 

Angular velocity 

Major axis half-length 

Axial ratios 

Bar angle 



Galactic bar 

0.62-109Mo 

60.0kms-^kpc-^ 

3.14kpc 

10:3.75:2.56 

20° 



(1) 





Spiral arms 




Number of spiral arms 


2 




Pitch angle 


15.5° 


(2) 


Amplitude 


3.36-10^Mokpc-3 




Angular velocity 


20.0kms-ikpc-i 




Scale-length 


2.6 kpc 


(3) 


Scale-height 


0.3 kpc 


(3) 


Scale-radius 


5.6 kpc 


(2) 



References. (1) iPichardo et aO (12004) : (2) iDrimmell (120001) : 
iJuricet all (120081) 



(3) 



and the angular velocity. For all paramet ers (listed in Table [2| ) 
we adopted the same values as reported by lPichardo et al.l (|2QQ4|) 
and refer the reader to their original paper for a justification of 
the various choices. 



2.3. Spiral arms 

The parameters of the spiral structure are much more uncertain 
than for the bar. Most observational studies suggest two-armed 
or four-arme d spiral patter n, or their super-position (see for ex- 
ain2le lValledl2005l:IValleell2QQ8l:lAllen et al.ll2QQ8l:lGrosbol et al.l 
l20lTh . 

We modeled the gravitation al potential perturbatio n induced 
by spiral arms as described by ICox & GomezI (|2Q02|) . We used 
a two-armed spiral pattern with their basic sinusoidal arms. The 
perturbation has the shape of a logarithmic spiral. The shape of 
the spiral arms model is defined by five parameters - radial scale 
length, scale height, pitch angle, amplitude of density perturba- 
tion (mass density of the arms in the Galactic plane at the solar 
radius), and angular velocity. The vertical structure of the arms' 
mass density has the form sech^(z/Zs), where Zs is a scale height. 
We chose the values of parameters consistent with a current pic- 
ture of the MW spiral structure - see Table |2] 

iDrimmel & Spergell (l200ll) present a three-dimensional 
model for the MW fitted to the CORE far-infrared and near- 
infrared emission. They found that two-armed structure domi- 
nates the near-infrared emission and the arm-interarm density 
contrast of their modeled stellar spiral arms has a value o f abou t 
1.8 at the solar radius (see Fig. 15 in Drimmel & Spe rgelll2QQll) . 
Our choice of parameters for Cox & G omez (2002) gives 2.0 
for the same ratio. The value is also consistent with the obser- 
vati ons of ex ternal galaxies with a weaker spiral structure (see, 
e.g.,|Blocketal. 2004). 

In addition to the parameters defining the shape of spiral 
arms, we needed to specify the initial orientation of the arms. 
We followed the result of the study of t he COBE/DIRBE K- 
band emission made by iDrimmell (l2QQQh . which is consistent 
with two-armed spiral model of the MW. 




Fig. 1. Maps of strength of non-axisymmetric components in 
the combined (bar and spiral arms) today's MW model (noted 
as MWl in Sect. |4l). Color scales map the ratio between ra- 
dial (left) and tangential (right) force resulting from the non- 
monopole part of the gravitational potential to the radial force 
due to the monopole part at the same radius. Maps are plotted 
for the Galactic plane. O symbol marks the position of the Sun. 



2.4. Strength of rotating patterns 

We estimated the contribution of the non-axisymmetric compo- 
nents by computing the force resulting from the bar or/and spi- 
ral arms' gravitational potentials in the disk plane, following the 
approach of Combes & Sanders (Il98ll) (see also Buta & Bloclj 
2001). Fig. [T] shows maps of the ratio between the radial, resp. 
the tangential, force due to the non-monopole part of the com- 
bined (bar and spiral arms) potential to the radial force due to 
the monopole part at the same radius. 

The strength of the bar is often characterized by the bar 
strength parameter, which is defined as the maximal value of the 
ratio between the bar's tangential force to the monopole radial 
force at the same radius. Our choice of bar model and parame- 
ters gives a bar strength parameter of 0.4. The model of spiral 
arms gives a value of 0.05 (using the tangential force from the 
spiral arms instead of the bar in the previous definition). 

2.5. Orbit calculation 

Once a model of gravitational potential is adopted, the cor- 
responding equations of motion can be constructed. We used 
Hamilton's equations, which give a set of first-order diff'erential 
equations that we integr ated using a Buli rsch-Stoer integrator 
with adaptive time-steps (JPress et al.lll992h . The relative change 
in the Jacobi constant for the models with a single rotating pat- 
tern is estimated to be on the order of 10"^^. 

To see and understand the diff'erent orbits caused by the bar 
or/and spiral arms we calculated orbits in four versions of the 
MW model - axisymmetric potential, potential with the bar but 
not spiral arms, with spiral arms and no bar and, finally, a com- 
bined model with both bar and spiral arms. 



3. Recent orbit from backward integration 

To investigate the shape and properties of the recent orbit of 
NGC 6791, we integrated its orbit backward for 1 Gyr, starting 
from the cluster's current position and velocity vectors. These 
initial vectors were obtained from the most recent observational 
estimates we found in the literature: equatorial celestial coor- 
dinates (a and S for J2000.0 equinox), distance to the Sun Jq, 
radial velocity v^, and PM. 
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Table 3. Input observational data used for the calculation of the 
initial conditions - equatorial coordinates a, S (J2000), radial ve- 
locity Vr, PM jda cosS, fis, and distance to the Sun J©; together 
with their references. 



Adopted value 



Reference 



a 


290.22083° 


(1) 


6 


37.77167° 


(1) 


Vr 


-47. 1+0.7 km s-i 


(2) 


fia COS 6 


-0.57+0.13 mas yr-i 


(2) 


1^5 


-2.45+0. 12 mas yr-i 


(2) 


d^ 


4.01±0.14kpc 


(3) 









References. (1) WEB DA; (2) iBedin et alJ ll2006il (see also references 
there-in); (3) lBrogaard et all (1201 Ih . 



3.1. Initial conditions 

As mentioned above, NGC 6791 has been the subject of intense 
study in the past. However, we rely in this study on the most 
recent values for the aforementioned parameters, which we list 
in Table [3] together with their source. So that they can be used 
as initial conditions for orbit integration, the input quantities re- 
ported in Table [3] need to be transformed into the appropriate 
coordinate system. 

We integrated the orbits in a cartesian galocentric right- 
handed coordinates system: the x axis points from the GC out- 
ward, in the direction of the Sun; the y axis points in anti- 
direction of the LSR motion; and the z axis toward the North 
Galactic Pole. The position of the Sun is given by vector 
(x,y,z)Q = (8.4, 0, 0)kpc in this coordinate system. The trans- 
formation of the observational data to the cartesian coordinate 
system centered on the Sun was per formed following step by 
step the iJohnson & SoderblomI (Il987l) algorithm updated to the 
International Celestial Reference System (the coordinates of 
the North Galactic Pole a = 192.85948°, 6 = 27.12825°). 
Adopting the solar motion with respect to the LSR (v© lsr) ^^m 
ISchonrich et aP (l201Qh : f ?7. V. W)r, = (11.1, 12.24, 7.25)kms-i 
(right-handed system, with U in the direction to the GC and V in 
the Galactic rotation direction) - the velocity vector with respect 
to LSR was obtained. Finally, using the Sun galocentric dis- 
tance of 8.4 kpc and the LSR rotation velocity of 243 km s"^ (see 
Sect. 12.1b , the position and velocity vectors were transformed 
into the coordinate system of integrations (denoted as x, y, z, Vx, 

The input observational data come naturally with some as- 
sociated uncertainties, which translate into uncertainties in the 
initial condition vectors. To take the observational uncertain- 
ties into account, we followed Dinescu et al. (1999) (see also 
IWu et al.ll2()()9l:IVande Putte et al.ll2Qia) . and calculated a set of 
1000 initial condition realizations. The initial velocity and posi- 
tion vectors were extracted from normally distributed values of 
radial velocity, PM, and distance to the Sun. The dispersions of 
a given distribution was taken as the observational uncertainty 
of the input data from Table [3] We integrated an orbit for each 
initial condition vector backward for 1 Gyr. The initial positions 
and velocities (x, y, z, Vx, Vy, v^) derived from the mean observa- 
tional input data are given in Table IH 



3.2. Orbits and orbital parameters 

Galactic orbits are routinely characterized by a set of so-called 
orbital parameters: the peri-galacticon Rp, the apo-galacticon7?a, 
the vertical amplitude |z|max defined as the maximal distance 



from the Galactic plane and, finally, the eccentricity e defined 
as (7?a - Rp)/(Ra -\-Rp). Depending on the type of orbit, values of 
orbital parameters can diff'er for individual revolutions. 

For irregular orbits orbital parameters change with time as 
the orbit evolves. This is particularly the case if the orbit is near 
a resonance for a single pattern or near a resonance overlap for 
bar-h spiral arms. Since chaotic behavior of orbits is expected 
in these regions, it is impossible to recover/predict an orbit on 
longer timescales. Therefore, to see the current orbital history, 
we calculated the orbital parameters for the most recent revo- 
lution (defined by azimuthal change of 27i) for each set of ini- 
tial conditions. The distributions of these orbital parameters are 
shown in Fig. [3] 

To highlight and quantify the influence of the various non- 
axisymmetric components, we integrated orbits in all four MW 
models described in Sect.[2l The orbits (backward integration for 
time of 1 Gyr) in axisymmetric and combined models given by 
the mean values of observational input data (see Table |4|for cor- 
responding initial conditions) are plotted in Fig. [2]- projection 
into the Galactic and meridional plane on the left and right, re- 
spectively. 

The orbits in models including only bar or spiral arms difl'er 
only slightly (around 10%) from the results of the axisymmetric 
model (see distributions of orbital parameters in Fig. |3). With 
increasing model complexity there is a trend of having larger 
apo-galacticon and lower peri-galacticon, which, in turn, yields 
a higher eccentricity. This is never higher than 0.33, however, 
which seems to make highly eccentric orbits for NGC 6791 im- 
plausible. 

According to these calculations, NGC 6791 is not expected 
to move much outside the solar circle, at odds with previous sug- 
gestions (see next Sect. l3.3l) . The orbit is presently located close 
to its peri-galacticon and this remains inside the solar circle on 
average, but never gets very close to the bulge. However, a sce- 
nario in which NGC 6791 formed close to the bulge and then 
migrated outward deserves a closer look, which we shall present 
in Sect. H 



3.3. Comparison with previous orbit integration 

To our kn owledge, the o rbit of NGC 6791 was stud ied three 
times -by ICarraro et all (l2006l) . iBedin et al.l (l2006h . and by 
IWu et al.' ('2009h. 

[Carraro et al] (l2006l) used c omplex potential derived from an 
N-body gas-dynamical model (iFuxll 19971 1 19991) . They obtained 
an eccentric orbit (e = 0.59) extending beyond the solar circle. 
However, the initial conditions were derived from absolute PM 
and radial velocity values available at that moment and admit- 
tedly less precise than the ones we are using here, and therefore 
it cannot be compared directly with the results presented here. 

iB edin et al.' ('2006') used a purely axisymmetric 
lAllen & Santillan (1991) model with the most up-to-date 
PM and radial velocity (the same values as we use here). 



However, reviewing the details of the initial conditions set-up, 
we identified a problem (mis-definition of the x axis orientation), 
which prevents a proper comparison with our results. 

■Wu et al.. (■2009.) calculated the orbit using the 
I Allen & SantillanI (1 19911) and two o ther axisymmetric po- 
tential models. As ICarraro et al.l (l2006h . they derived the initial 
conditions from less prec ise data than we use here, but still their 
orbital parameters for the lAllen & SantillanI (1 1991 1) potential are 
within error-bars consistent with results given in this study. 
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Fig. 2. Recent portion of Galactic orbit of 
NGC 6791 -orbit projections into the Galactic 
and meridional plane are plotted in the left and 
right column, respectively. Orbit in axisym- 
metric potential is plotted with the blue line. 
The red line shows the combined model (bar 
and spiral arms). The initial conditions (see 
Table lU are given by mean observational input 
data (Table[3l), the black dot marks the initial 
position. Integration was backward for time of 
IGyr. 
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Fig. 3. Distributions of cluster's recent orbital parameters. The distributions for peri-galacticon Rp, apo-galacticon/?a, vertical height 
of orbit Izlmax, and eccentricity e are plotted in columns. Different colors show results for different IMW models: axisymmetric 
model -light blue, model with bar -black, model with spiral arms -yellow, and red for the bar and spiral arms combined model. 



Table 4. Initial conditions for orbit calculations derived from the mean values of observational input data (Table [3]), together with 
standard deviations of their distributions. 



[kpc] 



[kpc] 



Z cr^ 

[kpc] 



[kms-1] 



[kms-^] 



7.05 0.04 -3.70 0.13 0.76 0.03 -39.7 2.7 -197.7 1.3 



Notes. See Sect. l3.1| for the description of the used coordinate system. 



[kms-1] 



-12.1 



Table 5. Parameters of initial condition distributions for the for- 
ward orbits integration. We used cylindrical galocentric coordi- 
nates - galocentric radius R, azimuth 0, vertical coordinate z and 
corresponding velocities vr, v^, v^. The tangential velocity v^p 
was calculated as the circular velocity Vc{R) in the axisymmetric 
potential as described in Sect. 12.11 at a given radius and then a 
random velocity with a dispersion of 6 km s"^ was added. 



Coordinate 


Distribution Center Dispersion, width 


R 

z 

Vr 


uniform 4 . kpc 1 . kpc 
0-2;7r with equidistant distribution 
normal Okpc 0.05 kpc 
normal Okms~^ lOkms"^ 
normal vdR) 6kms~^ 
normal Okms"^ 4kms~^ 



Table 6. Parameters for IMW models used for the forward inte- 
grations. 



Parameter 




MWl 


MW2 


Bar angular velocity 


[kms-^kpc-^ 


60.0 


45.0 


Bulge mass Mc2 


[lO^^Mo] 


0.62 





Bar mass 


[lO^^Mo] 


0.98 


1.6 


Major axis 


[kpc] 


3.14 


5.0 


spiral arms' angular velocity 


[kms-^kpc-i] 


20 


30 


spiral arms' amplitude 


[lO^Mokpc-3] 


3.36 


5.14 


spiral arms' pitch angle 


[deg] 


15.5 


20.0 



Notes. Both models have the same form and are characterized by the 
same parameters as the model used for forward integrations described 
in Sect. [2 (denoted as model MWl here). In model MW2, listed param- 
eters are changed (the rest of model parameters are kept the same as for 
MWl model, see Table[2l). 



In any case, the present study supersedes any previous one, 
in the sense that statistical analysis is performed for the orbital 
parameters adopting different Galactic model versions. 



4. Forward orbit integration and radial migration 

The inner disk, close to the bulge, is a high-density region 
wher e star for mation is very efficient and metal enrichment fast 
(Be nsby et al.i 2010). It is therefore conceivable to imagine that 
NGC 6791 could have formed there and then diffused outward. 
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Fig. 4. Time evolution of the apo-galacticon R^ against peri-galacticon Rp plane for the MWl model. A mean probability of finding 
an orbit within a given bin of /?a and R^ (where the mean is calculated over the indicated 1 Gyr time interval) is mapped. The bin size 
is 0.1x0.1 kpc. The probability for every bin is calculated as a sum of orbits with the bin's ^a and Rp values, where the contribution 
of each orbit is divided by the total number of revolutions per orbit within the time interval (so the total contribution of each orbit 
to each map is equal to 1). Finally, the probability is normalized to the total number of orbits 10^ and converted to %. Each column 
shows orbits for models including diff'erent components - axisymmetric model, model with bar, model with spiral arms, and finally 
bar and spiral arms combined model. Orbits are mapped within time intervals of 1 Gyr (0-1, 3-4, and 7-8 Gyr) as indicated in the 
left upper comer of each plot. The NGC 6791 's values of 7?a and Rp with their three sigma error-bars are plotted with a red cros s in 
each plot. The four lines in each of the plots show the lines of constant eccentricity for values 0.0 (circular orbits), 0.25, 0.5, and 
0.75, respectively, with decreasing slope. 
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Fig. 5. Time evolution of the apo-galacticon 7?a against peri-galacticon Rp plane for the MW2 model, see Fig. |4] for a detailed 
description of the plot. 



To estimate the probability and eflftciency of radial migration we 
carried out a forward orbits integration. More specifically, we in- 
tegrated forward in time a number of orbits with an initial loca- 
tion close to the GC. By comparing their orbital properties with 
the recent NGC 6791 orbital parameters obtained in Sect. |3] we 
attempted to estimate the probability that radial migration moved 
the NGC 6791 orbit outward from the inner Galactic regions. 



We followed 10"^ initially close-to-circular orbits distributed 
between 3 and 5 kpc from the GC. Initial positions of orbits were 
equidistantly separated in azimuth. Then, the radius was gener- 
ated randomly for each azimuth, from a uniform distribution be- 
tween 3-5 kpc (see Table \5\ for a detailed description of initial 
condition distributions in phase- space). 

To quantify the relative influence of difl'erent non- 
axisymmetric components, we followed the orbits in all four fla- 
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Fig. 6. Time evolution of the apo-galacticon 7?a against vertical amplitude |z|max plane for the MW2 model. Similar to Fig.lH results 
for the models including different components are plotted in columns; time increases with rows. 
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Fig. 7. Time evolution of the peri-galacticon Rp against vertical amplitude |z|max plane for the MW2 model. See Fig.|4]for a detailed 
description. 



vors of the Galactic potential model described in Sect. |2] (ax- 
isymmetric, barred, spiral, combined). When included, the am- 
plitude of the non-axisymmetric components grows simultane- 
ously in time for 0.4 Gyr from to its maximum value with the 
time-dependency described by Dehnen (2000, given by his equa- 
tion (4)). The gradual growth of non-axisymmetric perturbations 
is introduced to allow a progressive adaptation of initially nearly 
circular orbits. The growth of a real bar and real spiral arms is 
certainly much more complex and not simultaneous, but model- 
ing these time dependencies is beyond the scope of our model. 
We tested that our results are not significantly sensitive, apart 
from the first few hundreds of Myr, to the above gradual growth 
when compared to orbit integrations in which the bar and/or spi- 
ral arms were switched on abruptly to their full strength from the 
beginning. 



According to iMinchev & Famaeyl (l2010l) 3 Gyr should be 
sufficient for disk radial mixing. On the other hand, the age of 
the NGC 6791 cluster is estimated to be ~8 Gyr as mentioned in 
Sect. [TJ Therefore, we integrated the orbits forward in time for 
8 Gyr. Finally, the orbital parameters were calculated for indi- 
vidual rotations of each orbit as in the Sect. 13.21 



The properties of the bar and spiral arms that we used for 
the integrations of the NGC679rs orbit (as Hsted in Table |2]) 
are based on recent observations of the MW (see Sect. |2]i. We 
indicate this present-day MW model as MWl. From N-body 
simulations and observations of external galaxies at various red- 
shifts, we can infer the evolution of bar and spiral arms with 
time. For example, owing to the expected increase of stellar ve- 
locity dispersion with time, the disk becomes more stable and 
thus its spiral structure weakens. On the other hand, gas inflow 
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Table 7. Radii (in kpc) of the most important resonances for the 
Galactic bar and spiral arms MW models - corotation Rqr, inner 
Lindblad resonances Rilr, and outer Lindblad resonances Rqlr- 



Resonance 

^ILR 
^OLR 



MWl 



MW2 



bar spiral arms bar spiral arms 



3.6 
1.4 
6.9 



12.3 
2.8 
20.4 



5.1 
1.7 
9.4 



8.1 

2.2 
13.9 



Notes. Values for both non-axisymmetric components and both MW 
models are given. See Tables |2] and [6] for angular velocities and other 
model characteristics. The axisymmetric model described in Sect. 12.11 
is used for the unperturbed potential. 



can cool down the disk and the structure can strengthen again. 
Central mass conce ntrations (e.g., gas inflow) can decrease 
bar a mplitudes (e.g., lAthanassoula et al.ll2QQ5l: iDebattista et al.l 
l20Q6l) and even destroy bars, reforming later from the newly 
formed stars, but with a decreased strength and size (e.g., 
lournaud & Combes 2002; Combes 2008). 

Therefore, although we do not know the MW bar and spiral 
structure history, we can safely assume that those were much 
stronger in the past. For this reason we also followed orbits in 
a potential including stronger and faster spiral arms, as well as 
a stronger bar. This model is denoted as MW2. In Table|6]we list 
the values of the parameters for which MWl and MW2 diff'er. 
We decided to show an example of a model that is deliberately 
tuned to favor the radial migration caused by stronger bar and 
spiral arms and also by the choice of the rotating speed of the 
spiral arms pattern, bringing its resonances closer to the GC. We 
note that the parameter space (bar and spiral arms shape, extent, 
strength, and pattern speed) within the framework of our test 
particle integrations in analytically expressed potentials is very 
large and its systematic exploration was not the goal of the paper. 

In detail, the bar semi-major axis of the MW2 model was 
chosen to be longer, of 5 kpc (but the axis ratio for other two 
axes was kept as for MWl shape of bar) and the bar itself was 
more massive- so that on a time scale of 0.4 Gyr it replaced the 
whole mass of the bulge. As a consequence, the strength parame- 
ter of this bar is 0.5. As for spiral arms, the spiral pattern was set 
to have a higher amplitude and higher pitch angle. The angular 
velocity of the bar in this case was decreased to 45 km s"^ kpc"^ . 
This brings the bar's corotation radius to 5.1 kpc and so the ra- 
tio of bar's corotation radius to its semi-major axis is consistent 
with the value of 0.9-1.3 found in most bar red galaxies (see 
paragraph 6.5.1 in iBinney & Tremainel 120081 and references 
therein). On the other hand, the spiral arms angular velocity is 
increased to 30kms"^ kpc"^. The resonances of spiral arms are 
therefore closer to the GC where the orbits are originally dis- 
tributed. See Table |7] for resonance radii in both models. 

Before starting to comment more on the results, we recall 
the reader that we aimed to search for similarities between the 
present-day orbital parameters of NGC 6791 and the orbital pa- 
rameters of orbits that started close to the GC and then migrated 
outward. If there were similarities, we could claim that the sce- 
nario in which NGC 6791 could form in the very inner disk is 
plausible also from a dynamical point of view, and not supported 
only by the observed chemical characteristics of the cluster. 

The outcome of these experiments is illustrated in Figs.|4]-|9l 
Figs.|4land[5]map the evolution of mean probability of orbit with 
a given combination of 7?a and Rp. For the direct comparison 
with the recent orbital parameters of NGC 6791 we also indicate 
its values with the red cross in the maps. Figs. [6] and |7] map the 







Fig. 8. One of the Galactic orbits and its orbital parameters evo- 
lutions calculated in the MW2 model including bar and spiral 
arms -upper left plot shows the orbit in the Galactic plane; up- 
per right plot shows the time-evolution of the vertical amplitude; 
the bottom row shows evolutions of eccentricity, and peri- and 
apo-galacticon in the left and right plot, respectively. For one of 
its revolutions, this orbit fits the 3cr intervals for peri-galacticon, 
apo-galacticon, and vertical amplitude of the NGC 6791 orbit, 
these moments are marked with circle symbols in plots on the 
right. 



mean probability in the vertical amplitude |z|max against the 7?a 
and Rp plane for the MW2 model, respectively. As mentioned 
above, model MW2 was chosen to obtain an enhanced migra- 
tion. Comparing Figs. |4] and [5l we can see that this is indeed the 
case, illustrated by the higher values of apo-galacticon seen in 
Fig.O Moreover, some orbits have rotations with similar orbital 
parameters as we derived for NGC 6791 (the position of the red 
cross in the 7?a against Rp plane is within the area populated by 
some of the forward integrated orbits). 

From Figs. [6] and |7] we see that only the MW2 combined 
model can produce orbits that also show a vertical amplitude 
similar to NGC 6791. Only very few orbits (169 out of 10"^) can 
at least once per the total integration time of 8 Gyr reproduce 
the 7?a and Rp values similar to NGC 6791, and only some of 
these orbits (13 out of the 169) meet all three orbital parameters 
together. To investigate these orbits in more detail, we display all 
three orbital parameters together in Fig. |9]for the time interval 
of 7-8 Gyr of the integration. 

The experiment with forward integrated orbits illustrated in 
Figs.|4]to[9]helps us to provide the following considerations: 

• only model MW2 that includes bar and spiral arms produces 
the apo-, peri-galacticons, and eventually also the vertical 
amplitude similar to the backward integration of the current 
orbit (see the maps for the combined model in the rightmost 
column of Figs. [5]-[7l). 

• Fig.[8]shows an example of an orbit that migrated close to the 
recent NGC 6791 one. Initially, the orbit is close to circular 
with the radius of 4 kpc, and later starts to oscillate between 
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Ra [kpc] 

Fig. 9. Apo-galacticon R^ against peri-galacticon Rp plane for 
orbits integrated in the MW2 model during the time interval of 
7-8 Gyr. The range of /?a as-well-as of Rp is limited to the 3cr 
intervals of orbital parameters values of NGC 6791. Each revo- 
lution with this 7?a and Rp is shown by a point. The color of each 
point maps the vertical amplitude Izlmax- Revolutions that match 
NGC 6791 for all three orbital parameters (within their 3cr inter- 
vals) are marked with squared boxes. 



very eccentric and close to circular again but has a radius of 
about 7 kpc. 

from our experiment, we estimate the probability for the 
MW2 model including bar and spiral arms to produce the 
present orbit of NGC 6791 restricted to 2D Galactic plane to 
be about 0.40% (during the total integration time of 8 Gyr). 
This probability is calculated in the same way as for maps in 
Figs. |4]-[7l as a sum of orbits reaching R^ and Rp within 3cr 
intervals of NGC 6791, where the contribution of each orbit 
is divided by the total number of revolutions per orbit within 
the time interval of 1 Gyr (so the total contribution of each 
orbit within the time interval is equal to 1). Finally, the prob- 
ability is normalized to the total number of orbits 10^ and 
we sum the probabilities for time intervals of 1 Gyr within 
the whole integration time of 8 Gyr. If the vertical amplitude 
is taken into account, the probability has an even lower value 
of 0.03%. 



5. Summary and conclusions 

We investigated a scenario in which the old, metal-rich open 
cluster NGC 6791 formed in the inner disk, close to the bulge, 
and then moved outward because of perturbations induced by 
the bar and spiral arms. This scenario might indeed explain 
its unique properties, since star formation close to the bulge is 
strong and metal enrichment fast. NGC679rs large mass and 
compactness can explain why it has survived to the present time. 
We conducted our investigation using a numerical tool that 
allowed us to integrate orbits in different Galactic potentials. 
We calculated the recent (1 Gyr back in time) Galactic orbit 
of NGC 6791 in various MW models - axisymmetric, including 
bar and/or spiral structure. The observational uncertainties were 
taken into account in Monte Carlo fashion, calculating a set of 
1000 orbits with their initial conditions distributed according to 
PM, radial velocity, and distance observational uncertainties. We 
found that orbits in the axisymmetric model and models includ- 
ing only one rotating pattern differ only slightly. The orbit in 
the combined bar and spiral arms model is more eccentric with 



eccentricity up to -0.3, and has a higher vertical amplitude of 
-Ikpc. 

To investigate whether these orbital parameters can be 
achieved as a consequence of migration processes from the inner 
disk, we also followed a set of forward integrations (for 8 Gyr) 
with initial positions closer to the GC (with galocentric radii 
of 3-5 kpc). To this aim, we considered two different realiza- 
tions of the MW- present-day model: the same we used for the 
backward integrations (MWl), and models with rotating patterns 
more supporting the radial migration process (MW2). 

The purpose of this analysis was not to reproduce the cluster 
dynamical history precisely, but to estimate how efficient the mi- 
gration induced by the MW bar and spiral structure could be and 
how high the probability of the orbit realization within the lim- 
itations of our model is. We indeed found that our MW2 model 
that incorporates both the bar and spiral arm perturbations , can 
produce orbits with apo-galacticon and peri-galacticon similar 
to the actual values of NGC 6791. However, the probability of 
this scenario, as quantified from our experiment, is low, approx- 
imately 0.40%. Moreover, our migration scenario struggles to 
reproduce the vertical amplitude of NGC 6791. This might be 
an artifact of the Flyn n et al.l (Il996b potential, which does not 
represent the vertical structure of the disk well over the stud- 
ied range of galocentric distances, as well as of uncertainties in 
the modeled vertical structure of the bar and spiral arms. We are 
aware that this could also imply a problem with the migration 
scenario itself, question that is unfortunately insolvable within 
the limitations of the current model. 

Our study is unable to provide evidence for the inner origin 
of NGC 6791 and its later displacement due to radial migration 
caused by the bar and the spiral arms resonance overlap. This 
could be viewed as an argument in support of the external (ac- 
creted dwarf satellite) origin of NGC 6791, in line with the UV 
upturn argument mentioned in Sect. \T\ On the other hand, our 
simple model gives a non-zero probability of migration from the 
inner galaxy, which would therefore deserve additional investi- 
gation with improved models of the MW and with a more sys- 
tematically covering of mode l's parameter space. 

The recen t resu lts by iMinchev & Famaeyl (l2010h and 
iMinchev et al.l (l201lh show that the addition of a second per- 
turbation creates stochastic regions throughout the galactic disk 
and thus allows for a large-scale migration without resorting to 
the effect of short-lived transients. We note that for transient 
spirals this effect would be amplified. We expect that the ad- 
dition of a secondary spiral den sity wave, as found in N-body 
simulations of barred systems (lOuillen et al.l 1201 Oi Minchev 
et al., in preparation ), that may possibly exist in the MW 
(Naoz & Shavivl 120071) . will make the migration process more 
efficient. Therefore, the strengths of bar and spirals we used in 
this study may be the upper limits on what is required to push 
NGC 6791 outward. 

Future studies should strive to incorporate more realistic dy- 
namical models, e.g., considering the time evolution of the bar 
and spiral structure parameters, the effect of a gaseous compo- 
nent initially present in the disk, and continuous gas accretion. 
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